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Abstract. We investigate the Gross-Pitaevskii equation which describes an atomic 
Bose-Einstein condensate confined in an optical lattice and driven by a spatiotemporal 
periodic laser field. It is demonstrated that the exact Floquet states appear when the 
external time-dependent potential is balanced by the nonlinear mean-field interaction. 
The balance region of parameters is divided into a phase-continuing region and a phase- 
jumping one. In the latter region, the Floquet states are spatiotemporal vortices of 
nontrivial phase structures and zero-density cores. Due to the velocity singularities of 
vortex cores and the blowing-up of perturbed solutions, the spatiotemporal vortices 
are unstable periodic states embedded in chaos. The stability and instability of these 
Floquet states are numerically explored by the time evolution of fidelity between the 
exact and numerical solutions. It is numerically illustrated that the stable Floquet 
states in the phase-continuing region could be prepared from the uniformly initial 
states by slow growth of the external potential. 
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1. Introduction 

The atomic Bose-Einstein condensates (BECs) in optical lattices have stimulated great 
interests in both the many-body quantum effects and mean-field dynamics |5]- 

[7j. It has been suggested for applications in quantum interference [H |2], quantum 
information processing [HI IHl [ID] and matter- wave manipulation [HI [T2]. Using the 
mean-field theory, some exact stationary solutions for the systems of one- dimensional 
(ID) optical lattices were obtained [131 El IE], and their stabilities are discussed 
[T6| [T7] . Some nonstationary states of BECs in time-dependent lattices have also been 
investigated [T8]-[2T]. 

The Floquet states, a kind of nonstationary states, have been extensively introduced 
to understand the dynamics of various driven systems [221 [23]. Fo^' a linear Schrodinger 
system of T-periodic Hamiltonian, H{f,t) = H{f,t + nT), the Floquet theorem 
allows one to write its states as \l/(r, t) = U{f,t)exp{—iEpt) with periodic function 
U{r,t) = U{r,t + nT) and quasienergy Ep. The Floquet analysis is analogous to the 
Bloch analysis in solid state physics in which the states of spatially periodic system are 
written in terms of Bloch states and quasimomenta. Recently, the Floquet analysis 
has been introduced to systems of condensed atoms. For driven systems of BECs 
in double-well potentials, applying a two-mode approximation to the Gross-Pitaevskii 
(CP) equation, the Floquet states have been used to analyze the coherent control of the 
population self-trapping [221 123] . For the driven ultracold Bose atoms in optical lattices, 
which obey driven Bose-Hubbard models in full quantum theory, the Floquet states are 
applied to investigate the dynamical superfluid-insulator transition [21]. Below, we will 
generalize the familiar Floquet states in driven linear Schrodinger systems to the ones 
in a driven nonlinear Schrodinger system without any approximation. 

Spatial vortices [25] are fundamental objects of spinning, often turbulent, flow (or 
any spiral motion) with closed streamlines. They widely exist in different fields including 
phase singularities in optics [26] and circulating particles in superfiuids and BECs [27] . 
Recently, the conception of vortex in spatial domain has been generalized to the one in 
spatiotemporal domain. It shows that the superposition of two phase modulated optical 
beams generates the train of spatiotemporal vortices, which are periodical in space or 
time p8]. We expect such spatiotemporal vortices can also appear in an atomic BEC 
under some certain conditions. 

The stability analysis can provide useful information for preparation, control and 
application of a particular state. Mathematically, the instability appears when an 
initially small deviations do not keep bounded. It has been demonstrated that a blowing- 
up solution [29] appears in a BEC system and is related to the BEC collapse and 
instability [30]. For a linear quantum system, the fidelity between the unperturbed and 
perturbed states, a quantum Loschmidt echo, has been successfully used to analyze the 
stabihty [311 (3^. In which, the fast decay of fidelity corresponds to the instability of 
a quantum evolution. In the following, we will extend such a fidelity analysis to the 
nonlinear quantum systems of BECs. 
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In this article, we show how to prepare the exact Floquet states of an atomic BEC 
trapped in an optical lattice and driven by a spatiotemporally periodic field, and also 
analyze their stabilities. Utilizing the balance condition between the external potentials 
and mean-field interaction, we obtain a kind of exact Floquet solutions for the non- 
integrable chaotic system. The balance region of parameters is divided into the phase- 
continuing region and phase-jumping one, which are associated with the stable periodic 
states and unstable spatiotemporal vortex states embedded in chaos respectively. The 
stability and instability are confirmed by the corresponding fidelities between the 
exact and the numerical solutions. Our results suggest a method for suppressing the 
instability and preparing stable non-stationary states of the condensates. The considered 
condensate stabilization and preparation could be experimentally realizable. 



2. Model and exact Floquet states 

We consider an atomic BEC with strongly transverse confinement, and so that it obeys 
a quasi-lD GP equation, 

iHt = -^^xx + [^idl^P + V{x, t)]^, (1) 

where m is the single-atom mass and gid denotes the quasi-lD interaction strength [33], 
and V{x,t) stands for the external potential V{x,t) = Vocos^(/cx) + f{x,t) with lattice 
potential Vocos^(fcx) of strength Vq and driving field f{x,t) to be determined, and the 
wave vector k. 

To obtain the exact solutions of model (1) with the balance technique [131 IE]) we 
use the balance condition 

9iM^ + Vix,t) = EF (2) 

with Ep being a constant, which means that the nonlinear mean-field interaction is 
balanced by the external potential. Subsequently, the model (1) is simplified as 

ihipt = -—ipxT + Epil) (3) 
2m 

whose complete solution which contains some arbitrary constants can be easily obtained 
[15j and the corresponding integration constants are determined by inserting the 
complete solution into Eq. (2). Obviously, Eqs. (2) and (3) supports the exact Floquet 
state, 



9id V 9id 

X exp (^-^Ept^ , (4) 



^{x,t) 



where ui is the driving frequency, Ep is the Floquet energy, ^jEp/gi^, y—Vo/gid and 
\/—EpVq are real quantities, and a denotes a constant with value being either 1 or —1. In 
the absence of a potential, V{x,t) = 0, the solution has an uniform density distribution 
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IV'P = IJ'/gid for = J jJi/gid e with ji = Ep. In the case of Vq 7^ and / = 0, the 



balance solution (4) contains the stationary state ^/'(x,t) = ay — Vo/5'id cos(fca;)e~*^* for 
the chemical potential fi = u = hk'^/2m and zero Floquet energy Ep = 0. In general, 
the balance condition (2) and Floquet solution (4) require the driving frequency u and 
recoil energy E^ of the optical lattice satisfy the relation 

huj = Er = — — , (5) 
2m 

and the driving field has the form 



f{x,t) = cos /ex cos tut, Vi = 2ay —EpVo, (6) 

where the driving strength Vi is determined by the lattice depth Vq and the Floquet 
energy Ep. Given Eq. (5) and the form of potential, the second of Eq. (6) becomes 
the balance condition corresponding to Eq. (2). This condition indicates the relation 
between the driving strength Vi and frequency u for the requirement of the balance. 



The real values of ^ Ep/ gi^ requires that the attractive {gid < 0) and repulsive (5^1^ > 0) 
condensates have negative and positive Floquet energies, respectively. The driving field 
in Eq. (6) is a time-dependent laser standing wave, which can be formed from the 
linear superposition of two counter-propagating travelling waves [I9]. It is clear that 
the classical system governed by the total potential V{x, t) = Vq cos^ kx + Vi cos kx cos ut 
is chaotic and non-integrable [3ll [351 [36] , the corresponding GP system (1) is also 
non-integrable chaotic one even for the stationary state case 

Starting from the non- stationary states (4), one can selectively prepare different 
stationary states of expected chemical potential = u via selecting the driving 
frequency u and adiabatically switching off the driving field, namely decreasing Vi 
to zero very slowly. Under the balance condition (6) this adiabatic operation makes 
Vi = and Ep = such that Eq. (4) becomes the simple stationary solution 
ip{x,t) = a^—Vo/gid cos(fcx)e~*^* for /i = and a = ±1, which still obeys the 
balance condition (2) and additional Eqs. (3) and (5). It has been demonstrated that 
such a stationary state which has zero points possesses instability or undetermined 
stability for gid > or gid < respectively [17|. Different from the stationary 
states, the exact Floquet state for a driven condensate has a nontrivial phase which 
cannot be simply decoupled into a spatially dependent part from the super flow and 
a temporally dependent part dominated by the chemical potential. In contrast, this 
nontrivial phase includes a spatiotemporally dependent part from the super flow and a 
temporally dependent part dominated by the Floquet energy. This is resulted from the 
spatiotemporal dependence of the super flows in driven condensates. 

Defining the average number of atoms per well as = Jq \ilj\'^d{kx)/k = Tx{Ep — 
Vq/2) / {kgid), the Floquet energy can be expressed as Ep = kgidN/n + Vq/2. The real 



values of ^Ep/gn and y-Vo/5'id in Eq. (4) require Ep/g^ = kN/ix + VQ/{2gid) > 
and the lattice strength Vq satisfying, 

1^1 < K = 2kN\g,d\/n. (7) 
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Given the Ep and condition Vo/gid < 0. the balance condition in Eq. (6) gives 
Vi = 2ay/^WVo = a^-2{2kNgid/n + Vo)Vo. This formula and Eq. (7) confine the 
driving field strength 

\V^\ = ^2|yo|(K-|1^o|) < VJV2, (8) 

where the inequality is derived from the maximum condition d|yi|/d|Vo| = for a fixed 
Vc- Therefore, Vc and Vc/V2 are the supercritical values of the lattice depth \Vo\ and 
driving strength |Vi| for the balance solution. 

Due to cos{kx) = i(e**^^' + e^*'^'^), the exact Floquet solution (4) can be regarded as 
the coherent superposition between an atom standing wave formed by two counter- 



propagating plane waves and the background y EF/gidG:x.p{—iEFt/7i). Writing the 
macroscopic wave function as ^ = R{x,t) exp[i9{x,t)], the exact solution (4) implies 
the atomic-number density 

R' = W (9) 

= — (Ef — Vq cos^ kx — Vi cos kx cos ujt) 
9id ^ ' 

and the phase 

e = arctan[/m(^)/i?e(7/')] (10) 



I Vo I cos kx sin ujt Ep 

arctan '' j= —t. 

\Ep\ + J \Vo\ cos kx cos ut 



Clearly, the atomic density R^{x,t) has the same profile of the potential function and 
can be controlled by the external fields. The exact formula for the amplitude R and its 
finite form will raise the control precision. The non-integrability of the system means 
that its solution cannot contain the integration constants determined by the initial and 
boundary conditions. Therefore, the general solution cannot be given exactly. Of course, 
one may write an exact solution in terms of a superposition of an ensemble of complete 
orthogonal sets including Fourier modes, which could contain infinite terms. 

3. Spatiotemporal vortices 

It is well known that the vortices are associated with the phase singularities. The phase 
accumulation along a closed circle around the vortex core and the phase jump along 
a line through the vortex core are determined by the vortex charge. According to the 
relation between the phase and the velocity, v = ^V9, the velocity divergence appears 
when there is phase jump. In our GP system (1), the velocity field is defined as = 
with 

kVi 

^x = 7, ^^7 — -rsinkxsinut (11) 

2gidR^{x,t) 

denoting the first derivative (gradient) of the phase (10). The corresponding fiow density 
J = R^Ox = sin /cx sin u;t with amplitude being proportional to the driving strength 
Vi describes the Floquet oscillations of the system. Given the density distribution (9), 
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we know that the macroscopic wave function has zero-density nodes {Xzd, tzd) satisfying 

R'^iXzd, tzd 



0, i.e. 



cos(A;x 



zd 



Vi cos (ujtzd) ± JV^ [cos^ (utzd) - 1] 



2Vo 



(12) 



Because cos{kxzd) are real quantities, one can get 



cos{ujtzd) = ±1, cos{kx 



zd) 



2Vo 



COs{ujtz 



from the above formula. Thus the spatiotemporal coordinates 
density nodes read 

vr 



■^zd 



tr. 



n- 



X 



10 

1 

k 



± arccos 



-Vi cosnvr 



2Z7r 



2K) 

[xzd, tzd) for the zero- 

(13) 
(14) 



with non- negative integers n and I. The phase gradient (11) becomes infinite at those 



zero-density nodes when sin (kx^i^ ^ 0. The infinite phase gradient 6x{x^i,tn 



means 



that the phase is a step function across the zero-density nodes (x^;,t„). Using the 
relation between velocity and phase gradient, we know that the velocity singularity 
appears at these points. Due to the nonlinear mean-field interactions otherwhere, 
Eqs. (11) and (9) imply that the velocity singularity at the density nodes with 
Ep = V{xzd,tzd) is an effect of nonlinear resonance 



The conditions (8) and sin (kx^i^ ^ which is equivalent to cos^ (jcx 
give the following parametric region, 

\Vi\ < 2\Vo\ and VJ3 < \Vo\ < K, 



< 1 



(15) 



in which phase jumps may occur. In this parametric region, applying the I'Hopital rule 
to the phase gradient (11), one can find 

k 



lim 9x 



t) 



±— lim 

2 t^-tn 



ta.n{kx^i) cot{ujt) 



oo. 



However, the time derivative of phase, 

lim et{x^i,t) = uj/2 - Ep/h 



t^tn 



(16) 



(17) 



keeps finite at {x^i^tn)- This means that the phase jumps occur along the spatial 
direction but not along the temporal direction. Thus, the whole balance region is 
divided into two parts with or without phase jumps. The boundary between two parts 

- K/3, 



has \Vi\ = 2J\EfVo\ = 2\Vo\, \Vo 



\E, 



i^nl^ tn) 



0, siyi{kx^i) = and 



finite 9x{x^i,tn). The latter indicates that the exact Floquet solutions at the boundary 
are not phase-jump solutions. In the phase-jump region, the driving field strength \Vi\ 
is smaller than two times of the lattice strength 2|Vo|. 

Similarly, applying the I'Hopital rule to the following quantities, one can easily 
obtain 



lim Rx 



t) cxi?(x±,t„) =0, 
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Figure 1. The spatiotemporal evolutions of the exact Floquet states: density 
evolutions (a) and (b), phase evolutions (c) and (d), and velocity evolutions (e) and 
(f). Left column [(a), (c) and (e)] corresponds to parameters Ep/gid — 3fc, Vo/gid = 
—0.3k, Vi/gid = 1.8974fc in phase-continuing region \Vi\ > 2|Vo|- Right column [(b), 
(d) and (f)] is associated with Ep/gid — 0.5k, Vn/gid = —2k, Vi/gid — 2k in phase- 
jumping region \Vi \ < 2|Vb|. The red dots in (f) are the points where phase jumps and 
velocity divergence occur. The atomic number-density is normalized in units of k, 
and the phase contribution from the factor exp(—j-Ept) is eliminated. 



lim Rt{x^i,t) = ViuJi/JVo, 



lim i?xxla;„;, 

£ — >tn 



lim 9^^[x^t, 



^ t)(x R-^ 
± t)^R-\x^i,t 



- oo, 
oo. 



With these formulas, one has 



and 



lim iljt{x^i,t) oc Rt{x^i,tn) = constant, 



lim ^^{x^i,t) (X 9^{x^i,tn)R{x^i,tn) = constant. 



Since ip, ipt and i/jx are bounded, the phase-jump solution ip is a bounded solution rather 
than a blowing-up one |29j . 

In Fig. 1, we show the spatiotemporal evolutions of the exact Floquet states for 
different parameters. The first, second and third rows show density, phase and velocity 
evolutions, respectively. The left and right columns have the parameters out and in 
region of phase jumps, respectively. For simplicity, without loss of generality, we rescale 
the density in units of k and eliminate the phase from the factor expl—jEpt). It clearly 
shows that the phase-jump solution (right column) has nodes of zero density at which vr 
phase jumps occur along the spatial direction and the corresponding velocities approach 
to infinite. 
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The nontrivial phase structure around the singular points, that is, the circulation 
integral §{9xdx + 9tdt) = § dO = 2mT {n = 0, 1,2---) along closed spatiotemporal 
trajectories enclosing a singular point, reminds us the well-known Onsager-Feynman 
quantization condition for planar vortices [381 [39], since Eq. (16) implies that the integral 
around the spatiotemporal singular points is not equal to zero, and Eq. (10) allows the 
transformation 6 6 ± nn at any spatiotemporal point. Therefore, mathematically, 
the circulation around the singular points turns out to be quantized as a consequence 
of the form of field 6{x,t). This means that the nonzero quantized integrals along the 
closed trajectories around the zero-density nodes indicate the existence of the (1 + 1)-D 
spatiotemporal vortex. Due to the analogues between the Floquet and Bloch analysis, 
the above Floquet states for an atomic BEC in ID optical lattice with a spatiotemporal 
driving field are similar to the nonlinear Bloch modes in 2D periodic potentials |40j . 
where the phase-jump solutions for the relatively weak driving strengthes correspond 
to the vortex solitons [H]. This peculiar type of vortices are called as spatiotemporal 
vortices [28]. Then the zero-density nodes (x^;,t„) are vortex cores. From the phase 
distribution in the right column of Fig. 1, we find the phase accumulation along a circle 
around the vortex cores (x^;,t„) is ±27r . This means that the vortex charge is ±1 and 
so that the vortices at a pair of (x^;,t„) and (x~;,t„) are a vortex-antivortex pair. 

4. Stability analysis and state preparation 

In this section, we analyze the stability of the exact Floquet solution (4) for the nonlinear 
quantum system (1). In the sense of Lyapunov's stability, the instability entails that 
the initially small deviations grow in time without upper limit. Within the linear 
stability analysis [171 [IB] ? the deviations are governed by the linearized equations for 
the unperturbed solution. If the linearized equations have any unbounded solution in 
time evolution, the corresponding unperturbed state is unstable. Only if all perturbed 
solutions are bounded, the unperturbed solution is stable. Usually, it is difficult to 
explore the stability of a nonstationary solution. However, it is relatively easy to show 
its instability. We will demonstrate that whether the nonlinear resonance from the 
velocity divergence can cause the dynamical instability. 

The perturbations to the exact Floquet solution (4) can be taken in different forms. 
We consider the one of them, which has been used for stability analysis of stationary 
state solutions [17], in form of 

11)= [R{x,t) + eipi{x,t)]exp[ie{x,t) - iEpt/h], (18) 

where the small parameter e obeys \e\ <^ 1 and the perturbation correction ^/'i(x,t) = 
(j){x,t) + iip{x,t) including a real part (j) and an imaginary part ip. According to the 
linear stability analysis, the Lyapunov stability of solution requires that the unperturbed 
and perturbed solutions satisfy the spatially boundary conditions and all first-order 
corrections including = \J\4'\'^ + Iv^P for arbitrarily boundary conditions are 

bounded in the time evolution, and the Lyapunov instability is associated with the 
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unboundedness of any perturbed solution. Substituting the perturbed solution (18) 
into the original equation (1), we get the linearized equations 

h(Pt = Li^ - 50, (19) 
hipt= ~iL3(P + Sip), (20) 

where the operators Lj and 5* satisfyl^^ 



Li = 

^ 2m 



^2 n2 

9l{x,t) +jgidR\x,t) 



+ V{x, t) + hOtix, t), {for J = 1 and 3), (21) 



2 



S =^ 

2m 



2e^{x,t)^ + e,,{x,t)]. (22) 

Using these operators and the relation ip = Rexp{i6), we can rewrite the unperturbed 
system (1) as LiR = and hRt + SR = 0. 

From the linearized equations (19)-(22), we know that the bounded solutions (p and 
ip must be zero at the zero-density nodes {x^i,tn) which are singular points of 9^ and 
9xx- It is quite difficult to solve the perturbed equations (19) and (20) for and ip. 
However, we can show the instability from the singularity of first derivatives 0t and 
ipt. Mathematically, the linearized equations (19) and (20) are a couple of partially 
differential equations with general solutions containing some arbitrary functions, which 
are associated with various initial perturbations. For the solutions of spatiotemporal 
vortices, the linearized equations (19) and (20) are singularly differential equations with 
infinite 6x and 6xx at the zero-density nodes (a;^j,t„). Perhaps one may find a set of 
special solutions for the linearized equations (19) and (20) whose first time derivatives (pt 
and ipt are bounded. Nevertheless, since (t){x^i,tn) and (p{x^i,tn) cannot keep vanished 
for the arbitrary functions included in the general solutions (f){x,t) and ip{x,t), we 
cannot guarantee the boundedness of the first time derivatives in the singular Eqs. 
(19) and (20). Hence, under some initial perturbations the singularity of ^^^(x^;,t„) 
and 9xx{Xni,tn) will result in the divergence of (pt^x'^i.tn) and Lpt{x'^i,tn). That is, 
the linearized equations (19) and (20) have blowing-up solutions |29] in the phase- 
jump region which supports trains of spatiotemporal vortex-antivortex pairs. The 
blowing-up of the perturbed solutions could break down the integrability of linearized 
equations (19) and (20) and induce the jumps in 0(x^;,t„) and (p{x^i,tn)- Moreover, 
the jumping heights A(p{x^i,tn) and Aip{x^i^tn) may be uncontrollably large, because 
of the non-integrability of the linearized equations (19) and (20) and the singularities of 
9l{x,t)(f)(x,t) and 9l{x,t)(p{x,t) at the vortex cores (x^;,t„). Therefore, the blowing-up 
of perturbed solutions relates to the instabihty of the system. 

To confirm the above theoretical analysis for instability, we numerically integrate 
the GP equation (1) with the well-developed operator-splitting method. In our 
numerical simulation, we input the initial conditions as the initial values of the 
corresponding exact Floquet states and choose the same parameters used in Fig. 1 and 
the periodic boundary conditions ip{xmax,t) = 4'{—Xmax,t) of Xmax = 4. To show the 
difference between the exact Floquet solution and the numerical one, we calculate the 
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time evolution of fidelity F(t) between the exact solution ipexi^, t) without perturbations 
and the numerical one ipnum{x,t) with initial white noise, 



1 



i^numix^'t)i^ex{x,t)dx 



Nc{t) 

where the normalization function is defined as 



2 



(23) 



Nc{t) = j lljl^^{x,t)llJnumix,t)dx 

X / ^*eAx,t)i^exix,t) dx, (24) 



which ensures the fidelity satisfying < F{t) < 1. Here, the functions ip*^{x,t) and 
ip*^^{x,t) stand for the complex conjugates of functions ipex{x,t) and ipnum(yX,t). Due 
to the GP system (1) conserves the total number of atoms, we have Nc{t) = Nc{0). 
The fidelity describes the similarity between the two solutions, F = 1 means the two 
solutions are the same and F = means no similarity between them. So the fast decay 
of fidelity indicates that the initially small deviation from the unperturbed state grows 
fast. This is a signature for the instability of a quantum evolution [311 [32] . 



2 3 4 5 1234S6 




4 8 1Z 4 S 1Z 



Figure 2. The numerical evolutions for the same cases shown in Fig. 1. The left [(a) 
and (c)] and right [(b) and (d)] columns have the same parameters with the left and 
right columns of Fig. 1, respectively. The above [(a) and (b)] and below [(c) and (d)] 
rows show the density and fidelity evolutions, respectively. For the non- vortex state 
(left column) the density evolution is regular and periodic, and the fidelity perfectly 
keeps close to one. But for the solution of spatiotemporal vortices (right column) after 
a short period of time the density evolution becomes chaotic and the fidelity shows 
fast decay, which is a signature of instability. 

In Fig. 2, we show the numerical simulations correspond to the exact evolutions 
shown in Fig. 1. For the non- vortex solutions (left column), the density evolution from 
numerical simulation is regular and periodic, and the fidelity between the numerical 
and exact solutions perfectly keeps close to one. This means that the two solutions are 
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almost the same and the exact Floquet state maintains its stability under the numerical 
perturbation. For the solution of spatiotemporal vortices, the density evolution becomes 
chaotic after a short period of time, and the corresponding fidelity quickly decays to 
zero after a short period of time and then chaotically oscillates around a small number 
close to zero. This means the spatiotemporal vortices being the periodic state embedded 
in chaos and indicates the system losing its stability under the numerical perturbation. 
In order to avoid such an instability, we need adjust the driving field strength |Vi| and 
lattice depth Vq to satisfy \Vi\ > 2|Vo|. Due to the instability is related to the nonlinear 
resonance, we call such an instability as the nonlinear resonance instability. 




Figure 3. State transition and preparation via slow varying processes, (a) Transition 
from a stable Floquet state to a uniform state by slowly ramping down the external 
potential V{x,t) = A{t)[cos'^ kx + {Vi/Vo)coskxcosLot] from A = Vq to A — and 
then keeping A as zero, (b) Preparation of a stable Floquet state from an uniform 
state of density B?{x) — n/gid with — Ep — Vb/2 and phase 9{x) ~ 9o by slowly 
ramping up the external potential V{x, t) ~ A{t)[cos'^ kx + (Vi/Vb) cos /cx cos wt] from 
A = to A — Vo and then keeping A as Vo- In which, the values of Vb, Vi, Ep, k and 
uj are as same as the ones for (a) of Fig. 1, and the linear ramping processes occur in 
the period from t = Q to t = Ibir /lo and then the amplitude A keeps unchanged after 
t = Xb-K /lj. 

The numerical and analytical results show that the condensate in the non- 
stationary states is stable in the parameter region > 2|Vo|, Vq < Vc/3. Noting 
the balance condition |Vi| = 2^\EpVq\ and the average number of atoms per well 
N = ndEpl + I Vo|/2)/(A;|(7i^|), for a fixed the Floquet energy Ep is controlled by the 
driving field strength Vi and lattice depth Vq. Therefore, one can control the oscillation 
amplitude and flow density of the stable Floquet states by adjusting the driving field 
strength and lattice depth. To precisely prepare the Floquet states analyzed above, one 
can use the similar techniques of adiabatic growth which has been successfully used to 
prepare stationary states for systems of time-independent potentials [13]. However, due 
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to our penitential is time-dependent even ignoring the time-dependence of the amphtude, 
the slow varying processes used here is non-adiabatic. Because of the conserved total 
particle number, the state ip = i?e~*^* of uniform density distribution R^{x) = fi/gid 
with /i = Ep — Vo/2 and trivial phase distribution 6{x) = 6q {6q is a constant determined 
by initial conditions) can be obtained by slowly ramping down the external potential 
V{x,t) = A{t)[cos'^ kx + {Vi/Vo)coskxcosu;t] from A — Vq to A — 0. Inversely, from 
the uniform state of density R'^{x) — {Ep — Vo/2)/fi'id and trivial phase, one can get 
the exact Floquet state (4) for the stability region \Vi\ > 2\Vo\, Vq < Vc/3 by slowly 
ramping up the external potential V{x,t) — A{t)[cos'^ kx + (T^i/V^) cos /cx cosa;^] from 
A ^ to A = Vq. 

In Fig. 3, we show our numerical simulation for the state transition and preparation 
via slow varying processes. For convenience, we choose the same values of Vq, Vi, Ep, 
k and uj for (a) of Fig. 1, which describes a stable Floquet state. In panel (a), whose 
external potential is linearly ramped down from A = Vq to A = 0, we show the density 
evolution of a condensate evolving from the stable Floquet state [which is shown in panel 
(a) of Fig. 1] to an uniform state. Accompanying with the decrease of 1^41, the density 
difference between different positions becomes smaller and smaller. After t — IBtt/cu, 
there is no external potential, the density R'^/k — {Ep — Vo/2)/{kgid)±0-l = 3.15±0.1. 
This means that the relative density difference A(i?^)/(i?^) is about 3%. The fidelity 
between the final state and the corresponding exact uniform state is above 98%. In 
panel (b), whose external potential is linearly ramped up from A = to A = Vq, we 
show the density evolution of a condensate evolving from the state of uniform density 
R^{x) = {Ep — Vo/2)/gid and trivial phase to a stable Floquet state. Accompanying 
with the growth of l^l, the density oscillation becomes more and more close to the one 
of the exact Floquet state. After t — IStt/o;, the fidehty between the prepared state 
and the exact Floquet state for the corresponding parameters is also above 98%. The 
numerical results may supply a useful benchmark for preparing and manipulating the 
stable Floquet states in experiments. 

5. Conclusions and discussions 

In conclusion, for an atomic BEC in a one-dimensional optical lattice, we have shown 
how to prepare the exactly nonlinear Floquet states in a balance parametric region 
via a spatiotemporally periodic driving field. The exact Floquet solution describes the 
interference between an atomic standing wave and an uniform background. The balance 
condition requires that the sum of internal and external potentials is equal to the Floquet 
energy which is proportional to the atomic number per well and the depth of the lattice 
potential. It is shown that the phase of the exact Floquet solution may be continuous or 
piecewise continuous, depending on the ratio between the driving intensity and lattice 
depth. The phase-jumping solutions with piecewise continuous phase include the trains 
of spatiotemporal vortex- anit vortex pairs which are embedded in chaos. Applying the 
linear stability analysis, we have analyzed the stability and instability of the exact 
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Floquet states. The instability is related to the blowing-up of the perturbed solutions 
and the nonlinear resonance. The stable periodic states and unstable phase-jumping 
states have been numerically illustrated by the corresponding fidelities between the 
analytically unperturbed solution and the numerically perturbed solution. Dividing 
the balance region of parameters into the stability and instability subregions, we can 
selectively prepare the stable periodic states and unstable spatiotemporal vortex states 
by adjusting the driving strength and lattice depth to fit the corresponding subregions, 
and prepare the different stationary states by adiabatically switching off the driving 
field. 

We expect the above exact Floquet states in nonlinear quantum systems would 
stimulate experimental interests in investigating and stabilizing the non-stationary and 
stationary condensates. It would be helpful to explore the chaotic dynamics in such 
non-integrable systems of spatiotemporal periodic potentials like V{x,t) [35l [M]- For 
the non-integrable system of driven condensates, the exact Floquet state is completely 
determined by the parameters of external field. Therefore, we can prepare and control 
the spatiotemporal structures described by the Floquet solutions via choosing and 
adjusting parameters in the balance region. For a small ratio between the driving 
intensity and lattice depth one could explore the instability from velocity singularity 
via observing the breakdown of the periodic structure in density distribution [l2l [19] . 
Utihzing the slow varying processes, one can selectively prepare the stable Floquet 
states. Particularly, the results points out a general route to experimentally stabilizing 
and preparing the non-stationary states of the condensate. 

Additionally, with the similarity between the Floquet and Bloch analysis, the well- 
developed techniques for manipulating the atomic condensates with laser and magnetic 
fields offer an opportunity to study the analogue and difference between Floquet states 
and Bloch modes including spatial and spatiotemporal vortices. 
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